Typeset with jpsj2.cls <ver.l.2> 



Letter 



o 
o 



On 



43 

o 



Effective Field Theory of Triangular-Lattice Three-Spin Interaction Model 
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We discuss an effective field theory of a triangular-lattice three-spin interaction model defined 
by the Z p variables. Based on the symmetry properties and the ideal-state graph concept, we 
show that the vector dual sine-Gordon model describes the long-distance properties for p > 5; 
we then compare its predictions with the previous argument. To provide the evidences, we 
numerically analyze the eigenvalue structure of the transfer matrix for p = 6, and we check 
the criticality with the central charge c = 2 of the intermediate phase and the quantization 
condition of the vector charges. 
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The XY model consisting of the inner products of two 
neighboring planar spins with the 7L p symmetry-breaking 
field has been extensively researched; it provides the ba- 
sic understanding of more complicated models and also of 
real materials. Especially, for the two-dimensional (2D) 
case, the effective theory for the long-distance behav- 
iors was established based on the 2D Coulomb gas (CG) 
picture and the renormalization-group (RG) concepts, 
thereby enabling its complete understanding. 1 In this pa- 
per, we investigate its extension, i.e., the three-spin in- 
teraction model (TSIM) that was introduced a long time 
ago by Alcaraz et al. 2 ' 3 Suppose that (k,l,m) denotes 
three sites at the corners of each elementary plaquette 
of the triangular lattice A (consisting of interpenetrating 
sublattices A a , Ab, and A c ), then the following reduced 
Hamiltonian expresses a class of TSIM: 



U = -ft J ^2 cos (<p k + ip i 
(k,l,m) 



(1) 



Angle ifk = 2Ttrik/p, nu £ [0,p — 1] defines the Z p vari- 
able. For p = 2, eq. (1) is the exactly solved Baxter- 
Wu model with three-spin-product interactions, but for 
larger p, it becomes puzzling in its expression in terms of 
the spin variables. However, for p > 5 on which we will 
concentrate, its effective theory possesses a remarkably 
simple structure and can provide a unified understanding 
in a wide area of researches, including those on the 2D 
melting phenomena. Therefore, our aim is to formulate it 
based on a recent approach and to confirm its predictions 
quantitatively via numerical calculations. 

We shall begin with the symmetry properties. Adding 
to the translations and space inversions, the model is 
invariant under the global spin rotations: ip^ i— > ifk + 
Sp= a b c SieA (2imp/p)fik,li if t ne sublattice dependent 
numbers satisfy a condition n a + + n c — (mod p) , 2, 3 
This symmetry operation — denoted as (n a ,nb,n c ) — can 
be generated from two of the three fundamental opera- 
tions with the following minimal spin rotations: 

£ a : (l,p-l,0), R b : (O.l.p-l), R c : (p- 1,0,1). (2) 
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Here, note that these commuting operations satisfy some 
relations, e.g., RP = RcR^R^ = 1 (identity operation). 
At a sufficiently low temperature, this symmetry can be 
broken due to the three-spin interaction, and one of the 
ordered states is stabilized. An obvious one is {fk} = 0, 
and the others can be obtained by successively applying 
eq. (2). Consequently, we can observe a p 2 -fold degener- 
acy corresponding to the order of the group. 2 ' 3 

Here, we shall mention that our strategy to formulate 
the effective field theory is based on the recent develop- 
ment by Kondev and Henley, 5, 6 where the critical ground 
states of various 2D classical systems have been treated 
in a unified way. So, to make it concrete, let us consider 
the structure of the so-called ideal-state graph X for the 
present problem. We expect that, like the flat states in 
the interface models, the above p 2 ordered states should 
be specified by the locking points of a certain kind of 
field variable; thus, the dimension and the structure of X 
determine the number of components and the compact- 
ification region of the fields, respectively. By definition, 
each node of X represents one of the ordered states, and 
each link exhibits the neighboring two nodes connected 
by minimal spin rotations. From the above-mentioned 
symmetry properties, there are two commuting genera- 
tors, so the graph is located in a 2D space. Further, due 
to relations such as R c RbR a = 1, the shortest loop of 
links should form the triangle. These require X to be de- 
fined on the triangular lattice. On one hand, for a given 
value of p, the periodicity conditions, e.g., R% = I, de- 
fine the unit cell of the so-called repeat lattice TZ. From 
these, we can obtain X. In Fig. 1, we give an example for 
p = 6, where 36 ideal states are specified by the sublat- 
tice dependent numbers n a , n^, and n c (e.g., "000" at 
the corners represents {fk} = 0). 

A two-component field theory is clarified to be relevant 
to the present model in its continuum limit; we shall 
write it as 3>(x) S M 2 (x is the 2D real-space position 
vector). It then satisfies the periodicity conditions $ = 
$ + 2Tre a , where {e Q } (a — 1, 2) denote the normalized 
(non-orthogonal) fundamental lattice vectors of 1Z (see 
Fig. 2). For the representations of the physical quantities 
in terms of it is also important to see the relationship 
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Fig. 1. The ideal-state graph X for p = 6. Triplets of numbers ex- 
hibit 36 ordered states. The contour plot of the locking potential 
and the transformation eq. (3) (inset) are also given. 




Fig. 2. The repeat lattice TZ and its reciprocal lattice TV. The 
fundermental lattice vectors are also given by arrows. 

between the discrete symmetry operations for spins and 
the field transformations, <I> i— > <I> + 5$>; especially, the 
correspondence of the minimal spin rotations eq. (2) to 
the following 5& (see the inset of Fig. 1): 

f 2tt p 2tt 2tt ? 2?r 

i? a : — ei, i? b : ei H e 2 , i? c : e 2 , (3) 

p p p p 

which retain relations, e.g., RcR^R^ = 1. Our next task 
is to derive the Lagrangian density by which the low- 
energy and long-distance behaviors can be captured. For 
this purpose, we shall first focus on the low-temperature 
region, at which the relevant potential may perturb the 
kinetic energy part representing the spatial fluctuation 
of <&(x). As we shall see, the latter is also responsible for 
the description of the intermediate critical phase. 2,3 

Due to the periodicity of <fr, the vector charge M is 
quantized to take the values on the reciprocal lattice 
of 1Z, 1Z* (Fig. 2). Then, the local densities relating to 
the spin degrees of freedom may be expanded to the 
Fourier series by using the vertex operators as Q($>) = 
Smstc* ^mc iM '*. While the inner-product form ensures 
a coordinate invariance in R 2 , we look into its expression 
by using their elements to get some familiarity. Writ- 



Author Name 

ing the reciprocal lattice vectors as {e a } (a = 1,2), the 
inner product between M = m a e a and <fr = V2<j) a e a 
becomes \/2m a (f) a , where the covariant and contravari- 
ant elements satisfy m a <G Z and V2(f> a = V2(f> a + 2ir, 
respectively (the summation convention is used and the 
factor y/2 is for later convenience), since {e a } and {e a } 
are dual (i.e., e a ■ e@ = 5%). Also, the norm of M, for 
example, is given by ||M|| = y/m a m a , where the metric 
tensor is defined by g a P = e a ■ e 13 and the contravariant 
elements by m a = g a Pmp, as usual. 7 

To obtain an explicit form of the locking potential, 
the following three issues need to be addressed: 5,6 (i) 
Since the potential is a part of the Lagrangian density, 
it should be invariant under eq. (3). (ii) It is sufficient 
to keep the most relevant terms in the RG sense. Since 
the dimension of a vertex operator is proportional to the 
squared norm of its vector charge (see below), it is suf- 
ficient to keep those with the shortest ones, (iii) The 
point group symmetry of the triangular lattice for I, 
which stems from the sublattice and the spin-rotational 
symmetries, should be taken into account. Consequently, 
these requirements can be realized as the expression 

Vm = ^E||M||= P a. eiM '*> where a * = Vd 11 and h p 
are the lattice constant of 1Z* and the coupling constant, 
respectively. The summation is over the following six vec- 
tors: ipe 1 , ±pe 2 , and ±p(e 1 +e 2 ). In Fig. 1, we give the 
contour plot of V(<&) for p = 6 and h 6 = 1/2; we ob- 
serve that the points with the minimum value form the 
triangular lattice and each point corresponds to the ideal 
state. 

While the locking point specifies the ordered state, the 
spatial fluctuation of the fields becomes important with 
the increase of the temperature. For p > 5, the critical in- 
termediate phase is expected, 2,3 and it must correspond 
to the roughing phase of an interface model. 8 Thus, we 
can safely introduce the free-boson Lagrangian density 
for the fluctuations: 

£„[*] = ^d^M", (4) 

where g a p = e a ■ e fj and </> a = gapcj) 13 . The Gaussian 
coupling K plays the role of inverse temperature to con- 
trol the stiffness of the interface. The summation over i 
(= x, y) specifying the Cartesian component of x in the 
basal 2D plain is also assumed. 
At this stage, £o + A with 

||M||=pa* 

{ijp.o oc hp and ao is the cutoff of A) can describe the 
lower temperature transition to the ordered phase. On 
the other hand, for the transition to the disordered phase, 
we should next consider the discontinuity of $ by an 
amount of 27rN (N g K), which becomes frequent with 
the increase of the temperature. This topological defect is 
created by the vertex operator e lN e> , where is the dual 
field of $ and is defined as Kdi<f> = eijdj® (cjj is the 
antisymmetric symbol). £o is invariant under the trans- 
formation, K <-» 1/K and $ <-> 0, so the stiffness of the 
interface described by becomes larger with the increase 
of the temperature. Therefore, the potential can become 
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relevant to bring about a unique flat state corresponding 
to the disordered phase. If we write = V29 a e a , then 
the covariant elements satisfy \[2Q a = \/28 a + 2ir, be- 
cause the vector charge N is in 1Z. To obtain the explicit 
form of the locking potential, we can repeat the same 
argument as before, and arrive at the following form: 

^] = ^E^ NS . (•) 

l|N|| = l 

which gives the extremum on the lattice points of 1Z*. 
Consequently, we see that the vector dual sine-Gordon 
field theory C = Co + C\ + C2 provides an effective 
description of TSIM. In the remaining part of this let- 
ter, we shall check our prediction analytically and nu- 
merically. Here, we summarize the scaling dimensions of 
the operators on Co- The vertex operator with the elec- 
tric and magnetic vector charges (M, N) is defined as 
Om,n = c 1 ( M #+N ' e> \ which has the dimension im,n = 
± (if ~ 1 1 1 M 1 1 2 + ^ 1 1 N 1 1 2 ) . This formula supports our treat- 
ment that the vertex operators with the shortest vector 
charges were kept in the effective theory. In addition, 
the bosonized expression of the spin degrees of freedom, 
Sk = c liPk , is also important. It only contains the elec- 
tric charges and should reproduce eq. (2) on applying 
cq. (3). We can see that the sublattice dependent ex- 
pression, (S a ,Sb,S c ) = (O e i + e2,o,0_ e i,o,0_ e 2 jO ), sat- 
isfies the requirement and reproduces the ordered-state 
spin configurations on X. Furthermore, these exhibit 
that the condition to give nonvanishing spin correla- 
tions 2,3 is reduced to that of the vector charge neutral- 
ity. 9 For instance, since the two-spin (three-spin) prod- 
uct 5 p (x)5*(y) [S p (x)S a (y)S T (z)] becomes neutral for 
p = g (p 7^ a 7^ t 7^ p), the average with respect to Co 
takes a nonzero value only between the same sublattice 
(among different sublattices). 

Now, we are in a position to discuss the transitions. 
Although the free part Co is perturbed by £1.2 (dimen- 
sions are x Vp0 — 2p 2 /3K and x yol = K/2), these are 
both irrelevant for p 2 /3 > K > 4. To obtain the lower 
and higher temperature transition points (say 7l,h), we 
should resort to the numerical calculations, but these 
are considered to correspond to the above terminuses. 
Thus, let us first focus on the region near Tjj where C\ 
is irrelevant. Instead of performing the perturbative RG 
calculation [olq — > a e dl ) 7 we shall see that our problem 
can be related to the works by Halpcrin, Nelson, and 
Young, 10-13 where the defect-mediated 2D meltings were 
studied based on the Kosterlitz-Thouless (KT) RG argu- 
ment 9 (the so-called KTHNY theory). In fact, Alcaraz et 
al., making use of the vector CG representation with the 
long-range interaction, pointed out the relevance to 2D 
melting via the dissociation of dislocation pairs without 
the angular force. While our Lagrangian density Co + C2 
is in the local representation, we can also find its coun- 
terpart in ref. 13. Therefore, we confirm their assertion, 
and we shall discuss the transition in detail. As the Burg- 
ers vector characterizes the dislocation for the triangu- 
lar lattice, N in eq. (6) takes the six vectors, ±ei, ±e2, 
and ±(ei — e 2 ), so the three-point function of the local 
density £2 provides a nonvanishing value owing to the 
satisfaction of the neutrality condition. This results in 



a nonzero operator-product-expansion (OPE) coefficient 
among them; it yields the 2/o i(0 term in the /3-function 
for the coupling constant J/0,1 (0 (the fugacity for disloca- 
tions), while that for x(l) = 2 — x yo 1 (I) basically remains 
unchanged. 14 Consequently, we obtain the KT-like flow 
diagram (see Fig. 8 in rcf. 13). For T > Th, the correla- 
tion length behaves as £ oc exp[C/(T— Th) s ] characteriz- 
ing the disordered phase, while for T < Th the exponent 
of the spin correlation function varies as 77 = 2xs = 4/3-fC 
and ?7h = 1/3 at Th (without corrections). 

Second, we shall discuss the region near TL where C2 
is irrelevant. While having seen the p 2 ordered states by 
C\ , we can also find a similar type of transition in refs. 10 
and 11, where the transition from the "floating solid" to 
"commensurate solid" was discussed for the case that the 
substrate periodic potential is commensurate with the 
adsorbatc lattice. Indeed, they observed its similarity to 
the dislocation-mediated melting. For the present model, 
Alcaraz et al. obtained the RG equations based on the 
duality observed in the vector CG representation. 2 ' 3 On 
the other hand, we can also reproduce them and also the 
finite £ for T < TL and tjl = 4/p 2 at XL, etc, based on 
our local representation and OPE coefficients. 

To confirm the effective theory, we shall provide the 
numerical calculation results here. We consider the sys- 
tem on A with M (— > 00) rows of L (a multiple of 3) 
sites wrapped on the cylinder and define the transfer 
matrix connecting the next-nearest-neighbor rows. We 
denote its eigenvalues as X q (L) or their logarithms as 
E q {L) = — i In |A 9 (L)| (q specifies a level). Then, the con- 
formal invariance provides the expressions of the central 
charge c and the scaling dimension x q in the critical sys- 
tems as E g (L) ~ Lf — ttc/QLC, and AE q (L) ~ 2irx q /L( > . 
Here, E e (L), AE q (L) [= E q (L) - E g (L)], ( (= 2/^3), 
and / correspond to the ground-state energy, an excita- 
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Fig. 3. (a) 1/T vs. c (circles) and r\ (squares), (b) The M depen- 
dence of rr M ,o at 1/T SD ~1.252. Dotted curves exhibit ||M|| 2 /2-Fsr 
with K = 4\/3 along the path in the inset. 



4 J. Phys. Soc. Jpn. 



Letter. 



Author Name 



tion gap, the geometric factor, and a free energy density, 
respectively. 15-17 When performing the diagonalization 
calculations, we employ two of three spin rotations in 
cq. (2) [e.g., (-R a ,.Rb)] as well as the lattice translation 
and space inversion. This is because the matrix size can 
be reduced, and more importantly, discrete symmetries 
can specify lower-energy excitations. For instance, we can 
find the excitation level corresponding to the spin oper- 
ator 5 a (x) in the sector with indexes (e l27r / p , 1). 

In the following, we give the results for p = 6, which 
are obtained using the data up to L = 9. In Fig. 3(a), we 
provide the T dependence (in units of J/ks) of c where 
the region keeping c ~ 2 can be recognized. From this 
plot, we roughly estimate the KT-like transition temper- 
atures as 1/TL ~ 1.5 and 1/Th — 1.1 (see also refs. 2 and 
3). Next, we plot the T dependence of r\ along its second 
axis, which shows that r\ is the increasing function of T 
and takes close values to ?7l,h = 1/9, 1/3 around the es- 
timates of T Lj h (see dotted lines). In Fig. 3(b), we give 
the scaling dimension x M as a function of the electric 
vector charge M along the path depicted in the inset. As 
a representative for the critical region, we pick up the 
self-dual point by numerically observing the level cross- 
ing, x ye o — x yo y, we then perform the calculations at 
1/Tsd — 1-252. From the figure, we can verify that the 
dimension depends only on the norm of the electric vector 
charge. Further, despite the treated L being small — the 
size dependence is actually larger for those with longer 
vectors — the results (open circles) agree well with the 
theoretical formula (dotted curves) , which quantitatively 
supports our above argument. 

To summarize, based on the ideal-state graph ap- 
proach by Kondev and Henley, we have discussed the 
effective field theory for TSIM. Due to the symmetry 
properties, two kinds of vector charges, M and N, take 
the values of the reciprocal and repeat lattice points, 
and then provide the descriptions of the phase transi- 
tions. While the effective theory is given by the vector 
dual sine-Gordon model, we have seen its close relation- 
ship to the KTHNY theory, as pointed out by Alcaraz 
et al. We also performed numerical calculations to check 
some theoretical predictions based on the conformal in- 
variance. 

Finally, we make two remarks: (i) For the determina- 
tion of Tl,Hj we observed the deviation from c = 2. This 
could provide rough estimations, but a more efficient cri- 



terion is desired. For the c = 1 case, the universal ampli- 
tudes of logarithmic corrections are utilized for the anal- 
ysis of the excitation spectrum, which then offers level- 
crossing conditions for the KT-transition points. 18 ' 19 For 
this issue, C provides the basic framework to perform the 
one- loop calculations, (ii) Apart from the above restric- 
tion p > 5, it is also interesting to sec whether C can 
describe the transitions observed in TSIM with p < 5. 
In particular, for p = 2, the existence/nonexistence of 
the critical RG flow connecting Cq to the Z2 orbifold of 
a Gaussian model with c = 1 (the critical fixed point 
of the Baxter- Wu model), due to the competing relevant 
perturbations £1,2, may be an important issue. 20 We will 
address these issues in the future. 
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the Japan Society for the Promotion of Science. 
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